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ABSTRACT 

We present results of deep radio continuum observations of two galactic globu- 
lar clusters 47 Tucanae (47 Tuc) and oj Centauri (u Cen) with Australia Telescope 
Compact Array (ATCA). No statistically significant evidence for radio emission 
was found from the central region for the two clusters. However, both clusters 
show a 2.5o" detection near the center that may be confirmed by future deeper 
radio observations. The 3a upper limits of the radio observations is 20 and 40 
/xJy/beam for u Cen and 47 Tuc, respectively. By using the fundamental plane 
of accreting black holes which describes the relationship between radio luminos- 
ity, X-ray luminosity and black hole mass, we constrain the mass of a possible 
intermediate-mass black hole (IMBH) in the globualar clusters. We also compare 
our results with other globular clusters and discuss the existence of IMBHs in 
globular clusters. 

Subject headings: black hole physics — globular clusters: individual (47 Tucanae, 
ui Centauri) 



1. Introduction 

Globular clusters are well-known exotic objects factories because of their highly fre- 
quent dynamical interactions. It has been decades since the prediction w as made that 



intermediate-mass black holes (IMBHs) mig ht exist in globular clusters (e.g., IWyllerlll970 



Bahcall &: Ostrikerlll975uFrank &: Reeslll976l ). IMBHs provide a possible connection between 
stellar-mass black holes and supermassive black holes, with their mass ranging from 10 2 to 
10 5 M . There have been many conclusive evidences supporting the existences of stellar- mass 
black holes as in X-ray binaries and supermassive black holes as in active galactic nucleus 
(AGNs). However, the existence of IMBHs in globular clusters is still under debate. Although 
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theoretical works provide possibilities of their existences, observational evidences are limited. 
There are many formation channels of IMBHs. One is that they may be the products of core 
collapses of population III stars (e.g.. lFryer et al.ll200ll ). IMBH s could also form through run- 
away merger processes in some y oung dense globular clusters fjPortegies Zwart fc McMillan 
20021 ; iPortegies Zwart et al.ll2004l ). Moreover, it is possible that in globular clusters, mergers 
of compact objects, such as stellar-mass black holes through interactions would result in a 
central black hole with mass ~ 0.00 1 mass of its host cluster, which falls in the mass regime 



of IMBHs (IMiller &: Hamilton! |2002| ). As a result, searching for IMBHs in globular clusters 
is crucial for understanding formation histories of dense stellar clusters, and it may further 
help us study the stellar dynamics and construct cluster evolution scenarios. 

There is an established "M-cr" relation between the superma ssive black hole mass 



and velocity dispersion in the host galaxie s based on observations (ITremaine et al.l 12002 



Ferrarese fc Fordll2005l ; iGiiltekin et al.l 120091 1 . The correlation implies that the formation of 
a supermassive black hole is possibly related to the formation of a galaxy. If the M-cr rela- 
tion extends to the regime of dense stellar cluster systems, globular clusters may host central 
IMBHs with masses ranging from 10 2 to 10 3 M . One of the formation scenarios of AGNs 
is that they may be grown from 10 2 to 10 3 M formed from direct collapse of population 
III stars. If IMBHs do exist in globular clusters, the M-cr relation may link the formation 
mechanisms of IMBHs and AGNs. 

As f ar as observa t ions are concerned, detecting IMBHs is achievable through many 
methods. iFarrell et al.l ( 120091 ) identified an IMBH in the galaxy ES0243-49 by its extremely 
high X-ray luminosity and variability in the X-ray luminosity and spectra. For IMBHs in 
globular clusters, the most common way is measuring velocities dispersion of cluster stars 
and then applying dynamical a modelling technique to constrain the mass-to-light ratio of 
globular clusters. Another method is detecting possible radio emissio ns from the center o f 
globular clusters expected from the fundamental plane of black hole (jMerloni et al.ll2003l ). 
The dynamical modelling technique can constrain the central " dark" mass in globular cluster, 
while it fails to distinguish whether the "dark" mass is contributed from a single IMBH or 
a group of neutron stars, white dwarfs, or stellar-mass black holes. Moreover, during the 
past decade, because of the improvement of the instruments, detection of possible radio/X- 
ray emissions predicted by the black hole fundamental plane becomes possible. As a result, 
several deep radio continuum observations of globular clusters have been carried out to search 
for IMBHs in globular clusters. 

To date, discoveries of IMBHs candidates in some globular clusters have been claimed 
with different observational techniques (Table 1). The report ed dynamical mass of a pos - 
sible IMBH in Gl (in M31) is M BH = 1.8±0.5 x 10 4 M f lGebhardt et~aD I2OO2L l2005h . 
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while in M15 (in our Milky Way), M BH = 1-7 1.7 x 10 3 M flGerssen et al.l 120031 ). Nev- 



ertheless, the claims of Gl and M15 based on dynamical method have suffered subsequent 
challenges . Theoretical models sugge st that it is not necessary to invoke a IMBH to explain 



the data (IBaumgardt et al.ll2003aU bl). Furthermore, latest radio continuum observations on 
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2008; 


> discovered in Gl 


(Ulvestad et al. 


2007 


). Nev- 



ertheless, recent Ch andra observation of Gl favours the idea that the central object in Gl 



is an X-ray binary (IKong et al.ll2010l ). Evidences from different detection methods are still 



not so compelling. It is thus inconclusive if any IMBH exists in globular clusters. 

The two most massive and concentrated Galactic globular clusters, u Centauri (u Cen) 
and 47 Tucanae (47 Tuc), hav e been proposed to harbour a cent ral IMBH. By using Hubble 
Space Telescope (HST) dat a. Ivan der Marel fc Anderson! (120101 ) set a la upper limit of 1.2 
x 10 4 M Q in u Cen, while iMcLaughlin et al.l d2006h put an upper limit of 1000-1500 M 
on the mass of a possible IMBH in 47 Tuc. In addition, previous radio observations set 
a conservative upper l i mit of 1000 M^ and 2060 M for u Cen and 47 Tuc respectively 
(IMaccarone et al.ll2005l ; Ide Rijcke et al.ll2006l ). As a consequence, our goal is to set a tighter 
upper limit on the mass of IMBHs in u Cen and 47 Tuc to provide evidences for or against 
the existence of IMBHs in globular clusters. 

Here we show the results of our deep Australia Telescope Compact Array (ATCA) radio 
observations of u Cen and 47 Tuc. We organize our searches of IMBHs in u Cen and 47 Tuc 
in the following sections. In §2, we describe the accretion models we adopted to constrain 
the black hole mass. The radio observations, data analysis and results are described in §3. 
We then constrain the masses of the possible IMBHs in both clusters based on the radio 
emission upper limits also in §3. Finally, we discuss the existence of IMBHs in globular 
clusters in section 84. 



2. Detecting IMBHs with Radio Continuum Observations 



Recent studies on correlations between X-ray and radio properties of both X-ray binaries 
and AGNs suggest that the ratio of ra dio-to- X-ray power increases with black hole mass and 



decreases with accretion r ate (see e.g.. lFalcke et al.ll2004 



Gallo et al 



Merloni et al 



2003; Maccarone et al. 



J2003h is shown in equa- 



20031 ; iMerloni et al.l 120031 ) . The correlation found in 
tion (1), where -Fsghz is the radio flux at 5 GHz, Lx is the X-ray luminosity, Mbh is the 
black hole mass, and d is the cluster distance. 
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It is believed that the black holes in globular clusters are likely more massive than 
stellar-mass black holes, and are accreting at low fractions of the Eddington luminosity 
(or we would have detect the X-ray emission from the cluster center). Consequently, the 
correlation of the fundamental plane may suggest that detection of radio power would be an 
important indication of IMBHs existing in globular clusters. 

Un der the circumstan ces of uncertain X-ray luminosity based on observational con- 



straints, 



Maccaronej (120041 ) further assumes that the IMBH would accrete intra-cluster gas 



via a Bondi -Hoyle-Lyttlet on (BHL) proces s (IHovle k, Lvttletonlll94ll ; iBondi k, Hovlel Il944 



Bondilll952l ). As shown by lHo et al.l (120031 ) and lFender et al.l (120031 ). the accretion efficiency 
(e) could be as low as 0.1% of BHL accretion rate, assuming a radiative efficiency (rj) of 10%. 
The BHL accretion rate, Mbhl, is described in equation (2), where n is the gas abundance; 
Tis the gas temperature of the globular cluster and assuming a general value of 10 4 K for all 
globular clusters. Given the fact th at some globular c lusters contain a substantial amount 



of gas measured from pulsars (e.g., iFreire et al.ll200ll ). one can adopt the gas abundance 



derived from these globular clusters in equation (2) (see Table 1). 
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Based on the calculation in iMaccarone &: Servillatl (120081 . 120101 ) . the radiative efficiency 
r) is expressed as 0.5mc 4 /LEDi> This relation is valid for the black holes in the low hard 
state and the X-ray luminosity Lx oc m 2 . The m here is the mass accretion rate. Thus, Lx 
could be estimated by BHL accretion rate Mbhl, and expressed with the black hole mass 
Mbh, the gas abundance n, and the gas temperature T of the globular cluster, as shown 
in equation (3). If we combine all the above information, the black hole mass could be 
constrained simply by radio flux and cluster properties. 
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L x = Tj e Mbhl (ergss 1 ). 



(3) 



u Cen was observed previously with ATCA simultaneously at 4.8 and 8.6 GHz over 12 hr, 
while 47 Tuc was observed briefly (3.5 hr) with ATCA at 1.4 GHz. Both observations failed 
to detect any central radio sour ce and yielded a 3a rms noise level of 96 nJy and 225 //Jy for 
uj Cen and 47 Tuc respectively (IMaccarone et al.ll2005l ; Ide Rijcke et al.ll2006f ). With different 
combinations of BHL accretion fraction and gas abundance in globular clusters, the radio 
upper limits correspond to the most probable estimation for the up per limit of the IMBH 
mass in u C en and 47 Tuc is 2000 M and 1000 M , respectively (see IMaccarone fc Servillat 
20081 . 120101 . and references therein). The results could not exclude the possibility that there 
is no IMBH in u Cen or 47 Tuc, which motivated us to propose deeper radio observations 
on both targets. 



3. Radio Observations 

3.1. Data Analysis 

47 Tuc and uj Cen were observed by ATCA during 24 to 25 January 2010. The observa- 
tions were performed simultaneously at frequencies 5.5 GHz and 9 GHz with configuration 
6A (with the baseline ranging from 337 m to 6 km) and the upgraded Compact Array Broad- 
band Backend (CABB). The data was taken with the CFB lM-0.5k correlator configuration 
with 2 GHz bandwidth and 2048 channels, each with 1 MHz resolution. The primary calibra- 
tor used was 1934-638 for both globular clusters, while the phase calibrators were 2353-686 
and 1320-446 for 47 Tuc and u Cen, respectively. At the start of each observation, we ob- 
served 1934-638 for 10 min. The phase calibrator was observed every 15 min. We use MIRIAD 



flSault. Teuben. fc Wrightlll995l ) to analyse the data with standard processes. When loading 
the ATCA data into MIRIAD, we use atlod with options birdie, xycorr, rf iflag, and 
noauto, which represents flagging out the channels affected by the ATCA self-interference, 
correcting the phase difference between the X and Y channels, discarding any autocorrelation 
data, and automatically flagging out frequency bands that are known to be heavily affected 
by RFI. We then perform the standard data reduction steps, including bandpass, phase and 
amplitude calibr ations. When producin g the dirty maps, we use multi-frequency synthesis 



(MFS) method (ISault fc Conway! Il999l ) and natural weighting to suppress the noise. The 
effective on-source integration time of 47 Tuc is 11 and 9 hours, while the effective on-source 
integration time of u Cen is 18 and 16 hours for frequencies 5.5 GHz and 9 GHz, respectively. 
The field of view of ATCA with configuration 6A is ~ 10' and 5' for 5.5 GHz and 9 GHz, 
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respectively, which both cover the whole area within the half-mass radius of 47 Tuc (r^ = 
3'17) and co Cen (r/j = 5(00). The spatial resolution of the radio observations can reach ~ 
1" to 2". 



3.2. Results 

Besides background quasars, we did not detect any radio sources at or near the central 
region of co Cen and 47 Tuc with radio emissions higher than 3a. In naturally weighted maps, 
the rms noise level for co Cen is 7 (5.5 GHz) and 11 (9 GHz) /xJy/beam; while for 47 Tuc, 
it is 17 (5.5 GHz) and 20 (9 GHz) /iJy/beam. We combined the observations with different 
frequencies in order to obtain a lower rms noise level at frequency 6.8 GHz. As a result, the 
observations reached a rms noise level of 6.5 and 13.3 /xJy/beam, giving a 3a upper limit of 
20 and 40 /xJy/beam for co Cen and 47 Tuc, respectiv ely. The brightest mill isecond pulsar 
in 47 Tuc has a radio flux of ~ 370 /iJy at 1. 4 GHz (IMcConnell et al.ll2004l ). Assuming a 



spectral index a = -1.8 for pulsar spectrum (IKramer et al.lll998l ). it would have a flux of 



~ 21 /iJy at 6.8 GHz. Thus the millisecond pulsars in 47 Tuc are all under the detection 
thresho ld of our observation . We overlapped the positions of the X-ray sources detected in 
co Cen (jHaggard et al.ll2009i ). and there is no any radio source with their signal-to noise ratio 
higher than 3 matched to the 180 Chandra X-ray sources. 

Figure 1 shows the HST optical images overlaid with ATCA radio map contours of 47 
Tuc and co Cen. The radio position of 47 Tuc and co Cen has an estimated error of ~ 1.4" by 
0.6" and 1.1" by 0.7" respectively, derived by dividing the beam size by the signal-to-noise 
ratio. Although we did not find any statistically significant radio emission, there is a 2.5<r 
detection near the center of both globular clusters (~ 0.4" from 47 Tuc center and ~ 0.7" 
from co Cen center). We calculated the average source number within the error ellipse of the 
globular cluster center by computing the area ratios of the error ellipse to the field of the 
10" x 10" region and also the total number of the sources (with their signal-to-noise ratio 
larger than 2.5) inside the field of view of the 10" x 10" region. With the assumption of 
Poisson distribution, we obtained the probabilities of finding one or more radio sources inside 
the error circle, which is ~ 2% for both 47 Tuc and co Cen. That is to say, the probability 
of chance coincidence is quite low and the nearby source could really be the central radio 



source. 
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3.3. Constraints on IMBHs mass 



The 3cr upper limits of radio emissions could be applied to the mass constraints on 
IMBHs in globular clusters, using the fundamental plane of black hole activity described in 
section 2. We estimate the mo st probable mass and conserv ative mass of the IMBHs in uo 
Cen and 47 Tuc as described in iMaccarone fc Servillatl d2008h. with the accretion e fficiency e 
= 0.1% - 3% and the radiative efficiency described in IMaccarone &: Servillatl (120081 ) . The gas 
abunance in globular clusters cou ld be estimated e ither based on dispersion measurements 
of pulsars, like Ml 5 and 4 7 Tuc ( iFreire et al.ll200ll ). or based on the s tellar mass loss rate 
(Ipfahl & RappaportkoOlh . We adopt the gas abundance calculated from lPfahl fc Rappaportl 



( 1200 ll ) in our calculation for u Cen, while for 47 Tuc, we use the gas abundance from pulsar 
measurements as their conservative estimation. The resulting mass of IMBH is 1100-5200 
M and 520-4900 M Q for u Cen and 47 Tuc, respectively. We list the summary of recently 
observed globular clusters and our observation results in table 1. 



Discussion 



To date, several radio continuum observations have been performed on a few globular 
clusters in order to detect possible IMBHs in cluster center. There is a radio source detected 
in Gl. However, it is still under debate that whether the radio source in Gl is related to 



an IMBH or a LMXB ( IKong et al.ll2010l ). Except for Gl, no central radio source has been 
detected in the other globular clusters. Nevertheless, the non-detections did not rule out the 
existence of IMBHs in globular clusters. The radio flux upper limits are used to constrain 
the IMBHs masses, regard ing; certain assumptions of gas properties in globular clusters and 
accretion model. Based on lPfahl fc Rappaportl (120011 ). the gas abundance in globular clusters 
could vary from 0.1 to 1 H em" 3 . The accretion efficiency of the putative IMBH in Gl is just 
below 1% ( lUlvestad et al.ll2007l ). Moreover, according to the standard accretion disk theory, 
the IMBH i n glo bular clusters would be in an extremely low state, while as suggested in 
Cseh et al.l (120101 ) . both the accretion efficiency and radiative efficiency could vary in a wide 
range. The scatter in the black hole fundamental plane correlation would also contribute 
errors in IMBHs masses estimate. Combined all these uncertainties, the estimation of IMBH 
masses is actually mo del- dependent and under strong model assumptions. 

Despite of the large uncertainties in the estimation of IMBHs masses based on radio 
observations, there are discrepancies between the IMBHs masses estimated from accretion 
model constraints and dynamically inferred IMBHs masses, especially for the case of u Cen 
(see table 1). The measurement of IMBHs masses based on stellar dynamics is limited to 
the instrument ability, which could not achieve high enough spatial resolution, although it 
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provides a more accurate estimation on the IMBHs masses. We further compare the upper 
limits of IMBHs masses of some globular clusters to the M-cr relation of galaxies (figure 2). 
T he upper limit for the mass of the IMBH in 47 Tuc is consistent with the M-cr relation 
of iGiiltekin et al.l ( 120091 ). whi le for u Cen, the upper limit of the IMBH mass is within the 
region of the M-cr relation of iFerrarese fc Fordl (120051 ). As the true masses of IMBHs must 
be lower than the upper limits, u Cen, 47 Tuc and M15 do not show inconsistency of the 
M-cr relation. There is no reason that globular clusters should obey the same relation which 
is confirmed to exist in massive galaxies. However, the consistency may suggest similar 
formation mechanisms or contents between globular clusters and galaxies. If we could lower 
the upper limits of the IMBH masses in globular clusters, it may provide evidences of different 
formation mechanisms between globular clusters and galaxies. 

In summary, we estimate the 3cr upper limit for the mass of the central IMBH in uj 
Cen of 1100-5200 M , and for 47 Tuc, we estimate the 3a upper limit of 520-4900 M . 
The estimations strongly depend on the accretion model and the gas properties of globular 
clusters. We detect a 2.5cr radio emission near the center of both globular clusters. Future 
radio observations with a sensitivity improvement to ~ 5 fiJy may be able to confirm if the 
sources are real or further constrain the parameters of accretion model. 

This project is supported by the National Science Council of the Republic of China 
(Taiwan) through grant NSC 96-21 12-M-007-037-MY3 and NSC 99-21 12-M-007-004-MY3. 
The Australia Telescope is funded by the Commonwealth of Australia for operation as a 
National Facility managed by CSIRO. 
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Table 1: Recent Radio continuum observations on globular clusters 



Cluster 


Distance 


n H 


Tgas 


-Fr,5GHz 


AfBH.rad 


MBH,dyn 


name 


(kpc) 


(Hcm" 3 ) 


(Kelvin) 


G"Jy) 


(M Q ) 


(M Q ) 


u> Cen 


5.3 


0.044 


10 4 


20 


5200/1100 


12000 


47 Tuc 


4.5 


0.28/0.07 


10 4 


40 


4900/520 


1500 


NGC 6388 


10.0 


0.1 


10 4 


81 


1500/735 


5700 


NGC 2808 


9.5 


0.26 


10 4 


162 


8500/1800 


2700 


M15 


10.3 


0.42/0.2 


10 4 


25 


4900/700 


1000 


M62 


6.9 


0.41 


10 4 


36 


2900/600 


3000 


M80 


10.0 


0.21 


10 4 


36 


5300/1100 


1600 


NGC 6397 


2.7 


0.16 


10 4 


216 


4300/900 


50 


Gl 


780 


-1 


10 4 


28 


4500 


18000 



NOTES. — The radio flux -Fr^ghz is the 3a upp er limit except for Gl, which is a detec- 

or NGC 2808 from 



tion. The radio flux for NGC 6388 is taken from lCseh et al.l (120101 ) 
Maccarone fc Servillatl (hoosh: for M 15. M 60 and NG C 6266 from 



Bash et all (120081 ) 



for NGC 6397 from lde Riicke et all tood ): for Gl from luivestad et al.l d2007h. The dv 
namic al black hole mass - MRH.Hyn adopted for oj Ce n is from Ivan der Marel fc Anderson 



( 2010 ); for 47 Tuc from McLaughlin et al 



(120071 ): for M15, M62 and M80 come from 



(2006); for NGC 6388 from 



Lanzoni et al. 



Bash et all (120081 ): f or NGC 2808 and NGC 



6397 from M-a relation of lTremaine et al.1 j2002 l): for Gl from iGebhardt et all d2002h . 



The radio black hole mass M BH ,rad for NGC 6388 is from lCseh et~aD (hold ). 
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Fig. 1. — The HST R-band image overlaid with ATCA contours at frequency 6.8 GHz within 
the central 10" x 10" region of the clusters. The red contours are la, 2a and 3d noise levels, 
while the green contours are -la and -2a noise levels. The center of the cluster is denoted 
by the blue cross. The blue ellipse represents the 95% positional error ellipse of the globular 
cluster center. Left : 47 Tuc. Right : to Cen. 
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Fig. 2. — The M-a relat i on of galaxies derived by iFerrarese fc Fordl ( 120051 ) (black-dashed 
line) and iGiiltekin et all ( 120091 ) (blue-dotted line). The red data points represents the 3a 
upper limits on the masses of IMBHs in glob ular clusters, with accretion efficie ncy e = 3%. 
The velocit y dispersion for u Ce n come s from Ivan der Marel fc Anderson bud); for 47 Tuc 



comes from lGebhardt fc Fischer! (119951); f or Ml 5 comes from (iGerssen et al. 



6397 comes from Mevlan &; Mayor! ( 1991 ); for Gl come s from Mevlan et al. 



comes from lMerritt et all (]200ll ); for Ml 10 comes from lValluri et all ( 120051 ). 
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